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by Robert G. Ragsdale 
Lewis Research Center 

SUMMARY 

An analytical study was made of the isothermal mass and momentum transfer that 
occurs when a heavy, slow- moving gas is injected coaxially into a duct of light, fast- 
moving gas. The outer stream is composed of two regions: an outermost region moving 
at high velocity, and an annular sheath of intermediate velocity that is immediately ad- 
jacent to the central core, or jet, inner fluid. As a basis for comparison, the total mass 
flow rate of outer fluid is held constant; thus, as the buffer region velocity is decreased, 
the velocity of the remaining outer region is increased. The effects of the velocity and 
thickness of this buffer region on the mixing of the two fluids are investigated for both 
laminar and turbulent flow. 

Outer- to inner-fluid initial-velocity ratios of 10, 50, and 100 are considered. Some 
representative velocity and concentration fields are shown for laminar and turbulent flow. 
The effect of a momentum buffer region is presented in terms of the amount of inner 
stream fluid that is contained within a specified length of the outer duct. The diameter 
of the duct is taken to be four times that of the initial jet; duct lengths of 2 and 4 jet 
initial diameters are considered. 

In a gaseous-fuel nuclear rocket engine, stream mixing is undesirable because it 
tends to dilute the central core of fissionable gas with the surrounding hydrogen propel- 
lant, and thus leads to an increase in required reactor pressure in order to maintain 
criticality. The results of this study show that M optimum M values of both buffer velocity 
and thickness exist which give a maximum amount of inner fluid in a given length of duct. 
The optimum velocity is greater for turbulent flow than for laminar flow. The optimum 
buffer thickness was found to be relatively independent of whether the flow is laminar or 
turbulent. For the system studied, the optimum buffer thickness is about 1 inner stream 
radius. The presence of a buffer region increases the amount of inner gas contained in 
a duct length of 2 initial jet diameters by more than a factor of 2. 

*Presented at the AIAA Propulsion Joint Specialists Conference, Colorado Springs, 
Colorado, June 14-18, 1965. 



INTRODUCTION 


The interaction between two coaxially flowing fluids is a subject of considerable 
interest, both academic and practical. A description of mass, momentum, and energy 
transfer between dissimilar coaxial streams is required in order to understand and 
design combustion chambers, jet pumps and ejectors, and afterburners. Recent experi- 
mental and anlytical studies have been reported that directly involve such a flow pattern. 
In a supersonic combustor (refs. 1 to 4) a hydrogen jet is injected into a coflowing stream 
of oxidizer gas. In a proposed concept for a nuclear rocket engine (refs. 5 to 7) a low- 
velocity, fissionable gas is injected into a coaxially flowing, high-velocity stream of 
hydrogen propellant. In both situations, the goal of the studies is to understand, to pre- 
dict, and, ultimately, to influence the mixing rate of the two species. 

For the gaseous -fueled- reactor concept, the mixing is undesirable, because it tends 
to dilute the nuclear fuel with hydrogen and thus leads to an increase in required reactor 
pressure. Flow schemes that result in slower stream mixing are therefore of interest. 
Since the mixing occurs primarily because the outer stream velocity greatly exceeds 
that of the inner stream, it is possible, in principle, to delay the interaction of the two 
streams by separating them by means of a buffer region of fluid flowing at some inter- 
mediate velocity. This report describes an analytical study of the characteristics of 
such a flow field. 

In particular, the system to be considered is one in which a jet of low-velocity, 
heavy gas is separated from a surrounding, high-velocity, light gas by an annular sheath 
of intermediate-velocity, light gas. The gases are contained within a duct as illustrated 
in figure 1. The situation considered is one in which the total mass flow rate of the 
light gas is held constant as it might be in a rocket engine that is required to produce a 
certain amount of thrust. Thus, as the velocity of the buffer region is reduced, the 

velocity of the remaining outer stream 
must be increased to maintain a constant 
mass flow rate. 

The numerical solution of refer- 
ence 5 and its extension to turbulent flow 
(ref. 7) have been used to investigate the 
effects of a momentum buffer region on 
the coaxial mixing of dissimilar gases 
for isothermal flow. The primary param- 
eters are the buffer region thickness and 
velocity. Calculations were carried out 
for both laminar and turbulent flow. 

Figure i. - Calculation model. Three values of average outer- to inner- 
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stream velocity ratio were considered. The results are compared with the case where 
no buffer region is present. The results presented indicate that a momentum buffer 
region can be used to alter significantly the rate of mixing near the injection point of a 
coaxial jet. 


SYMBOLS 


b width of mixing region, ft 

c concentration of inner stream species, mole fraction 

D dimensional diffusion coefficient, ft /sec 

D dimensionless diffusion coefficient 

2 

D 1 1 molecular self-diffusion coefficient of species 1, ft /sec 

? 2 
Dj g binary diffusion coefficient of species 1 and 2, ft /sec 

K constant 

M molecular weight, lb/lb mole 

Re^ q initial Reynolds number of inner stream, 2r ^ qU^ Qp ^ q/^ 1 0 
r radial coordinate, ft 

r dimensionless radius, r/r ^ q 

Sc^ q initial Schmidt number of species 1, p ^ q/Pj q^I \ 
u axial velocity, ft/sec 

u dimensionless axial velocity, u/u^ q 

Vg molecular volume ratio 

v radial velocity, ft/sec 

v dimensionless radial velocity, v/u^ q 

z axial coordinate, ft 

z dimensionless axial position, z/r^ q 

/3 molecular weight parameter, (Mj/ft^) - 1 

9 

e eddy diffusivity, ft /sec 

e + dimensionless eddy viscosity, pe/p 

n K normalized containment factor 
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rj containment factor 

p local viscosity 

71 dimensionless viscosity 

p density, lb/ft^ 

\}/ stream function (see eqs. (5) and (6)) 

Subscripts: 
av average 

b buffer 

c l centerline 

max maximum 

min minimum 

t total 

w wall 

0 injection point (z = 0) 

1 inner stream species 

2 outer stream species 
Superscripts: 

— dimensionless quantity 
= average over entire outer fluid region 

ANALYTICAL METHOD 

The basic analytical method used was reported in reference 5 for laminar flow and 
is extended to turbulent flow in references 6 and 7. It will be briefly reiterated here, 
along with a more detailed discussion of its application to a turbulent three- region sys- 
tem. 


Laminar Flow 

The steady- state boundary-layer equations for isothermal, constant pressure, 
axisymmetric flow are 
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Continuity: 


— (pvr) + — (pur) = 0 
9r 9z 


Momentum: 


Diffusion: 


9u 


+ u 


9u 


9r 9z 


JL JL 

pr 9r 



v 3 £ + u = 

9r 9z r 9r 


rD l, 2 9c\ 
V P 3r / 


( 1 ) 


( 2 ) 


(3) 


It is assumed that there are no significant pressure gradients due to friction or momentum 
transfer. The dimensional equations (1), (2), and (3) are made nondimensional by nor- 
malizing all quantities to a dimensionally similar quantity in the inner stream at the in- 
jection point. For example, the dimensionless viscosity p is the ratio of local viscos- 
ity to that of pure inner stream fluid p/p^ q; the axial velocity is normalized to the jet 
injection velocity so that u = u/uj q; and the axial position and radius are normalized 

to the jet radius, so that z = z/r^ q and r = r/r^ q. The dimensionless diffusion co- 

efficient is given by D = (D^ 2^1 1^’ w ^* ere l the self-diffusion coefficient of 
the inner stream fluid. 

The dimensionless continuity equation is now 

— Fr v (/3c + 1)1 + — [ru(/3c + 1)1 = 0 (4) 

9r L J 9z L J 


where /3 a (M^/Mg) - 1. The system of equations is transformed from the r, z-plane to 
the ip, z-plane, where the stream function ip satisfies the continuity equation (eq. (4)): 

A. = _ (rv)(/3c + 1) (5) 

9z 


— = (ru)(/3c + 1) ( 6 ) 

9r 
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Introducing dimensionless quantities and the stream function along with the usual 
boundary-layer assumption that dxp/dz, « dip/ dr yields the following equation set: 



2 

0z 


r* dr* = / dhp* 

J 0 u(/3c + 1) 

^ + i_(ilur 2 (/3c + 1) iH.1 

Re 10 a^L 


(7) 

( 8 ) 


0c = 2 (/3c + 1) _l[du7 2 ( 9 ) 

3z Re^ q® c 1 o 


The boundary conditions are 


3 \js dxf/ 

= “ 2, 0 


z > = 0 


The local dimensionless viscosity is calculated from a mixing law utilized in refer- 
ence 5 that can be written in terms of dimensionless quantities as 


- = M 

03 + i)c + 

H 2 


( 10 ) 


The dimensionless diffusion coefficient is obtained by writing the Gilliland Equation 
(ref. 8) in the form 


D s 


D l, 2 _ 2 V2(/3 + 2) 

D ^M /3 ) 2 


( 11 ) 
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where Vg is the ratio of molecular volumes, as given in reference 8, of species 2 to 
species 1. 

Equations (7) to (11) are solved numerically by a computer program (ref. 5) to obtain 
radial profiles of dimensionless velocity and concentration at specified axial positions 
downstream from the injection point. In addition to constants relating to the numerical 
integration scheme, the following inputs are necessary and sufficient to define a case: 

(1) Physical properties: j8, JJ[ V 2 , SCj Q 

(2) Flow parameter: Re^ q 

(3) Initial velocity profile: u(0, r) 

(4) Initial concentration profile: c^T) 

For all the cases discussed in this report, slug flow velocity profiles are used, as 
illustrated in figure 1 (p. 2). For concentration, the initial profile is 

0 < r < 1 c = 1 


1<r - r max c = 0 

With a buffer region, the initial velocity profile is 

0 < r < 1 u = 1 


1 < r < r b u = u b 


iv < r < r u = u on 

b max 2 , 0 


and with no buffer region is 


0 < r < 1 u = 1 


! < r < r max u = u 2,0 


Turbulent Flow 

The laminar analysis is extended to include turbulent flow by adding turbulent trans- 
port coefficients to their molecular counterparts. Thus, the contribution of turbulence 
can be expressed in terms of a turbulent viscosity pe and a turbulent diffusivity e. 

It is also assumed that the eddy diffusivities for mass and momentum transfer are equal. 
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The total transport properties to be substituted in equations (2) and (3) in place of 
the laminar properties are written as 


P t = p + pe 

(12) 

D t = D + e 

(13) 


A similar procedure is employed in references 1 to 3, where it is further assumed that 
pe » p and e » D. 

Assuming that the same mixing law as for equation (10) can be used to obtain a 
dimensionless total viscosity and diffusion coefficient in terms of a dimensionless eddy 
viscosity e + results in 


^t = 


/3c + 1 

(/3 + l)c (1 - c) 

1 + €"*" (1 + € + )jU2 


(14) 



/3c + 1 

(/3 + l)c + 1 - c 


D + e+SCj 0 D + e + p 2 Sc i + 1) 


(15) 


where 


e + -£i 

M 

The functional dependence of e + remains to be specified in terms of spatial flow 
and/or property characteristics of the system. This has been accomplished in current 
coaxial mixing studies by utilizing Prandtl* s postulate that in a region of free turbulence 
the eddy diffusivity is proportional to the width of the mixing region and to the velocity 
difference across it: 


e = K^max - u min) < 16 ) 

This postulate and its application to single- component jets and wakes are discussed 
in reference 9. 

The extension of equation (16) to variable-density systems can be achieved in dif- 
ferent ways. It has been suggested that pe is constant with radius (refs. 1 and 2) and 
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that the velocity difference should be replaced with |P 2 U 2 - “ p cjt u cJ ( ref> Refer- 
ence 4 proposes that pe is constant with radius and that a combination of mass flux pu 

2 

and momentum flux pu be employed. In reference 10 it is assumed that pe/p is 
independent of radius, but dependent on axial position to some exponent. Reference 7 
concludes that pe/ p. can be taken as independent of radial and axial position and that the 
velocity difference should be expressed as (lug/u^ - l| )^ 2 . 

Although the various approaches have yielded unresolved differences, they all share 
some significant similarities. Each of these previous studies (1) report general agree- 
ment with some experimental data, (2) appear to be applicable for initial velocity ratios 
both above and below 1, (3) indicate that the radial variation of pe is not significant, 
and (4) conclude that density ratio must be included as a parameter when applying equa- 
tion (16) to variable -density systems. Each approximation has been shown to correlate 
data within certain ranges of flow and property parameters. This suggests that the dif- 
ference between laminar and turbulent flow is much greater than the differences between 
the various algebraic descriptions of turbulent viscosity. The numerical solution of 
reference 5 and the turbulent viscosity expression of reference 7 have been used for the 
calculations reported herein because no linearization assumptions are required and be- 
cause general agreement with average radial concentration data has been demonstrated 
(ref. 7) over the widest range of initial-velocity ratios (0. 83 < u 2 0. 97, 

1.25<u 2 <49). 

A recent study (ref. 11) has shown that these various correlations predict about the 
same turbulent viscosity at the jet origin when they are compared at the same velocity 
ratio. It was also shown that, for nearly equal stream velocities, turbulent mixing can 
be reduced by the presence of honeycomb inserts to remove initial turbulence from the 
two streams just prior to the jet origin. The turbulent viscosity correlation of refer- 
ence 7 used here included both the initial turbulence and that due to the shear between 
the two free streams. 

From reference 7, the ratio of turbulent to laminar viscosity can be written as 

£+ ‘7 =012 S)© <Rei ’ 0 ' 250,l " 2 ' 1|1/2 <17) 

For a two- region coaxial-flow situation, u ^ is simply the ratio of outer to inner 
initial fluid velocity, and the use of equation (17) is straightforward. When a momentum 
buffer is present, there are two velocities in the outer fluid, and it must be assumed 
that equation (17) still applies, with some kind of an average velocity used to get an 
"effective” Ug. It is not clear how to obtain one velocity to represent all of the outer 
fluid. Certainly this velocity should be some kind of average of the buffer region velocity 
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and the outer region velocity, but there are a number of possibilities, for example, a 
simple numerical average, a mass -flow- rate- weighted average, or a momentum- weighted 
average. The latter two procedures seem more realistic than a numerical average. For 
the most extreme case considered, the buffer region contained only 10 percent of the 
total mass flow of the outer fluid and 4 percent of the total momentum. Thus, either of 
the two weighted averaging methods produced a velocity that was essentially the same as 
that of the high-velocity outer region. The choice, then, was rather arbitrary; the Ug 
in equation (17) was evaluated with a momentum- weighted average velocity of the outer 
fluid. 


SCHEDULE OF CALCULATIONS 


The primary variables to be considered are the thickness of the buffer region, the 
velocity ratio in the buffer region, and the ratio of light-gas to heavy-gas initial velocity. 
Calculations were made for both laminar and turbulent flow. 

The first of each series of cases was for no buffer region present. For example, an 
initial-velocity ratio of 10 was selected for the first case. For the next case, a buffer 
region thickness of 0. 5 initial inner stream radius ( r^ = 1.5) was assigned, and the 
velocity ratio in this region was decreased until it was some specified fraction of that of 
the remaining outer fluid. A dimensionless channel radius r of 4 was assigned for 
all calculations, and as the buffer velocity was decreased, the velocity of the remaining 
outer fluid was increased to maintain a constant total-mass-flow rate of light gas. Thus, 
if the initial velocity ratio with no buffer region is Ug q, then 



For the illustrative case selected, Ug q = 10 and r^ = 1. 5; for a buffer velocity that is 
0. 1 of the remaining outer fluid velocity, 



10. 8 


and u b = 1. 08. 

Average outer- to inner-velocity ratios of 10, 50, and 100 were considered. For 
each of these ratios, buffer- to outer -velocity ratios of 0. 1, 0. 2, 0. 3, and 0. 4 were 
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investigated. For all combinations of these variables, buffer thicknesses of 0. 5 and 1 
initial jet radius were calculated, as well as the situation of no buffer region. The 
following combinations were calculated for both laminar and turbulent flow: 


Average outer- to inner-velocity ratio, Ug g 10, 50, 100 

Buffer- to outer-velocity ratio, u^/tig 0, 0. 1, 0. 2, 0. 3, 0.4 

Buffer region radius, r^ 1, 1. 5, 2. 0 


For all calculations, the following constant values were used: /3 = 100, 7^2 = 0. 2, 

V g = 0. 2, and Sc^ q = 1.0. As pointed out in references 5 and 10, the analytical results 
are relatively insensitive to the transport property ratios, so these choices are not crit- 
ical. Therefore, the results of this study should be reasonably general and not restricted 
by the particular properties chosen. Some of the calculations performed for this study 
were repeated for molecular Schmidt numbers of 0. 5 and 2, and the concentration and 
velocity profiles were virtually unaffected. Insensitivity of a similar analysis to Schmidt 
number was reported in reference 4. The laminar calculations are based on an inner 
stream Reynolds number of 200 and the turbulent cases for a Reynolds number of 20 000. 


RESULTS AND DISCUSSION 

The results of the calculations of three-region, two-component mixing will be pre- 
sented in two main categories. First, detailed characteristics of concentration and ve- 
locity fields will be discussed for both laminar and turbulent flow. Second, a more 
generalized presentation of the effects of a buffer region will be made in terms of the 
total mass of inner stream fluid present in a channel of specified length. 


Flow Field Characteristics 

The general effect of a faster moving outer stream is to accelerate and dilute the 
inner stream fluid. This results from mass diffusion of the inner fluid radially outward 
into regions of higher velocity and also from acceleration of the inner fluid by momentum 
transfer radially inward. Figure 2(a) shows the concentration field for turbulent coaxial 
flow with an initial velocity ratio of 50. If the initial velocity ratio were lower than 50, 
the mixing of the two streams would occur more slowly. A buffer region of lower ve- 
locity immediately adjacent to the inner stream presents such a situation, at least ini- 
tially. 

Figure 2(b) shows the effect of reducing the velocity in a buffer region 1 jet radius 
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Figure 3. - Centerline concentration. Turbulent flow; 
u 2 o = 50; buffer thickness, 1. 


Dimensionless 
axial velocity, 



(a) No buffer region. 


u 2 f o " 58# 



(b) Buffer thickness, 1; u^^ q = 0. 3. 

Figure 4. - Effect of buffer region on velocity field. Turbulent flow; u 2 0 = 50. 
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in thickness until it is 0. 1 that of the remaining outer stream. In order to maintain the 
same mass flow of outer fluid, or the same average velocity ratio of 50, the outer-fluid- 
velocity ratio between the buffer region and the channel = 4 is increased to 61 and 
the buffer- velocity ratio is then 6. 1. Figure 2(b) illustrates that the general effect of a 
momentum buffer is to reduce the mixing rate in the region near the jet. Farther down- 
stream, the mixing zone produced by the shear between the high velocity outer flow and 
the slow moving buffer region penetrates to the inner zone and results in increased ac- 
celeration of the inner fluid. 

Figure 2(c) shows the concentration field that results if the buffer- to outer-velocity 
ratio is increased from 0. 1 to 0. 3. Although the buffer velocity is higher, the mixing 
rate is less than for a buffer-velocity ratio of 0.1. This indicates that there is an opti- 
mum buffer velocity and that buffer velocities below this value do not afford sufficient 
momentum in the buffer to offset the increased momentum in the outer stream required 
to maintain a constant mass flow. A comparison of the concentration lines c = 0. 6 in 
figure 2 shows this effect. 

Figure 3 is a plot of centerline concentration as a function of axial position for the 
three cases first discussed. The centerline concentration is slightly higher for a buffer- 
velocity ratio of 0. 1 for only 0. 5 jet radius downstream; beyond 2 radii downstream it 
falls a little below that for no buffer region. Between axial distances of 1 and 4.5 jet 
radii, a buffer-velocity ratio of 0. 3 yields the highest centerline concentration. Fig- 
ure 3 also serves to illustrate that centerline variations are relatively insensitive to ef- 
fect of a buffer. 

The effect of a buffer region on the velocity field is quite pronounced. This is shown 



Figure 5. - Effect of buffer region on inner fluid average 
velocity. Turbulent flow; U 2 g = 50-, buffer thickness, 1. 



Figure 6. - Effect of buffer region on total mass of inner 
fluid in system of length z. Turbulent flow; U 2 g - 50; 
buffer thickness, 1. 
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in figure 4, again for turbulent flow and an average outer- to inner-velocity ratio of 50. 
Figure 4(a) is for no buffer region, and figure 4(b) is for a buffer thickness of 1 jet radius 
at a buffer- to outer-velocity ratio of 0. 3. The reduction in the rate at which the inner 
fluid is accelerated is clearly shown by a comparison of the constant-velocity lines for 
u = 5 (fig. 4). The constant -velocity lines for u = 49 show that farther downstream the 
higher velocity outer fluid begins to overcome the initial buffer effect. 

Figure 4 shows the entire velocity field regardless of species and does not, there- 
fore, clearly reflect what happens to the inner stream fluid. Figure 5 shows the axial 
variation of the radially averaged velocity of the inner fluid, still for turbulent flow and 
for a buffer thickness of 1 jet radius and an average initial-velocity ratio of 50. The 
effect of a buffer region is shown more clearly in figure 5 than by the centerline concen- 
trations shown in figure 3. A buffer- to outer-velocity ratio of 0. 3 reduces the rate of 
acceleration of the inner stream fluid. At a distance of 8 jet radii downstream, the 
presence of a buffer with a velocity ratio of 0. 1 results in a higher velocity inner fluid 
than the case with no buffer region. 

The influence of a buffer region on the average inner fluid velocity can be interpreted 
as an effect on the mean residence time of the inner fluid in a given length of channel. 
Thus, the effect of a buffer region can be expressed in terms of its effect on the total 
mass of inner fluid that is contained in a channel of length z, compared with the maxi- 
mum amount that could be present. This maximum would occur for slip flow and no 
radial diffusion. The ratio of the total amount of inner fluid present to that which would 
exist with no interaction between the two streams is then a sort of ” containment effi- 
ciency”, or containment factor. It represents the degree to which the inner stream 
fluid has been accelerated. It can vary from a maximum of 1. 0 down to l/ug q, which 
would represent instantaneous acceleration of the inner stream fluid. 

Figure 6 shows a plot of this containment factor rj c for various channel lengths 'z, 
turbulent flow, a buffer thickness of 1 jet radius, and ^ g= 50. Figure 6 shows the 
significant effect of a buffer region on the amount of inner stream fluid contained in a 
channel of length "z. For a system that is 4 jet radii in length, a buffer region with a 
buffer- to outer-velocity ratio of 0. 3 increases the amount of inner fluid present from 
18 percent without a buffer to 31 percent of the theoretical maximum. A buffer-velocity 
ratio of 0. 1 yields a containment factor of 0. 26. 

Figures 7 and 8 show turbulent concentration fields for initial-velocity ratios of 
10 and 100, respectively. Comparison of these figures with figure 2 (for Ug q = 50) 
indicates that the effect of a momentum buffer is more pronounced at higher initial- 
velocity ratios, quantatively shown in succeeding figures. 

For an initial-velocity ratio of 50, the ratio of turbulent to laminar viscosity is 830 
with no buffer region; this value was used in the calculation illustrated in figure 2(a). 
With a buffer region present, a momentum- averaged velocity of the outer fluid yields a 
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u 2( o = 12 - 2 * 



(b) Buffer thickness, 1 ; u^U 2 q = 0. 1. 


Figure 7. - Effect of buffer region on concentration field. Turbulent flow; 
a 2j o = 10. 
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u 2,0 * 100 "^ 



\\ , , . 

c o. 8-* \-.6 V4 V2 V 02 


(a) No buffer. 


u 2, 0 = 120 * 



(b) Buffer thickness, 1 ; u b /u 2 q = 0.2. 

Figure 8. - Effect of buffer region on concentration field. Turbulent 
flow; u 2 0 = 100. 
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(a) No buffer region. 




u 2,0 = 61 * 



(b) Buffer thickness, 1; u^/ u 2 0 a 0.1. 

Figure 9. - Effect of buffer region on concentration field. Laminar flow; U 2 0 = 50. 



value of 917 for pe/p for the case shown in figure 2(b). Figures 9(a) and (b) show the 
corresponding concentration fields for laminar flow (pe/p. = 0) for these two cases. A 
comparison of the concentration lines for c = 0. 2 in figures 9(a) and (b) shows that a 
momentum buffer also has a considerable effect for laminar flow, as for turbulent flow, 
even though the general nature of the flow fields for laminar and turbulent mixing is 
markedly different. 


Effect of Buffer Parameters on Containment Factors 

Containment factors for laminar and turbulent flow are shown in figure 10 for basic 
coaxial mixing, with no buffer region present, and for initial-velocity ratios of 10, 50, 
and 100. As before, the containment factor is the ratio of the total mass of inner fluid 
contained in a duct of length z to the amount that would be present with no interaction 
between the two streams. This containment factor serves as a measure of the degree of 
stream mixing. Higher initial-velocity ratios, longer ducts, and the presence of turbu- 
lence all lead to increased mixing and lower containment factors. 

The presence of a momentum buffer region reduces stream mixing and increases the 
containment factor. The effects of buffer- region thicknesses and velocities on the con- 
tainment factor are shown in figure 11 for laminar flow and in figure 12 for turbulent 
flow. Curves are shown for buffer-region thicknesses of 0. 5 to 1.0 jet radius and for 
initial-average-velocity ratios of 10, 50, and 100. 

The curves exhibit the same general effect of buffer velocity on containment factor 



Figure 10. - Containment factors. No buffer region. 
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(a) Axial position, ilr^ q = 4 ; aver- (b) Axial position, z/r^ q = 8; aver- 
age initial velocity ratio, u 2 g = 10. age initial velocity ratio, u 2 g = 10 * 



Buffer- to outer-velocity ratio, uj/u 2 q 


(e) Axial position, z/r^ g = 4 ; aver- (f) Axial position, z/r^ g = 8; aver- 
age initial velocity ratio u 2 g 55 100. age initial velocity ratio, u 2 g = 100. 

Figure 11. - Buffer containment factors. Laminar flow. 
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(a) Axial position, z/rj 0 - 4 ; aver- 
age initial velocity ratio, 0 - 10. 


Buffer to jet- 
radius ratio, 

r b 



(b) Axial position, z/rj 0 * 8 ; aver- 
age initial velocity ratio, u 2 g * 10. 



age initial velocity ratio, u 2 g = 50. 



(d) Axial position, zlr± 0 = 8; aver- 
age initial velocity ratio, u 2 0 = 50. 


.4 



0 .4 .8 1.2 0 _.4_ .8 1.2 

Buffer- to outer -velocity ratio, u b /u 2 q 


(e) Axial position, z/r 1 Q = 4 ; aver- (f) Axial position, z/r x 0 = 8; aver- 
age initial velocity ratio, u 2 g = 100- a 9 e initial velocity ratio, u 2 g * 100. 
Figure 12. - Buffer containment factors. Turbulent flow. 
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e 



(b) Buffer- to jet-radius ratio, r b = 2.0. 


Figure 13. - Normalized buffer containment factors for 
z,r i,o* 4 - 



Figure 14. - Op imum-buffer to primary-velocity ratios for z/rj g « 4. 
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Figure 15. - Optimum buffer thicknesses for 


z/r 


1,0 


■4. 


for laminar and turbulent flow. As the buffer- 
region velocity is decreased, the containment 
factor increases to a maximum value and then 
begins to decrease. A comparison of two cor- 
responding curves for z values of 4 and 8 shows 
that the buffer effect is reduced for longer ducts, 
but that the optimum buffer- to outer -velocity 
ratio is essentially the same, at least at the 
higher initial-velocity ratios. The remaining re- 
sults will be presented for a duct that is 4 jet 
radii in length (z = 4). 

It is difficult to assess the relative effects 
of a buffer region from the various curves of 
figures 11 and 12 because the absolute value of 
the containment factor is strongly affected by the 
initial velocity ratio and by the nature of the flow. It is therefore useful to normalize 
each curve to its value with no buffer region (u^/u^ q = 1.0). 

The variation of such a normalized containment factor 77^ with buffer- to outer- 
velocity ratio is shown in figure 13. The parameter 77^ is the factor by which the 
amount of inner fluid in the duct has been increased by a buffer region. For an average- 
velocity ratio of 100, a buffer that is 1 initial jet radius thick ( r^ = 2), and turbulent flow, 
the amount of inner fluid in the duct is increased by a factor of more than 2 by decreas- 
ing the buffer-velocity ratio ^ from 1 (no buffer) to 0. 2. The maximum buffer 

effect is obtained for a buffer-velocity ratio of approximately 0. 2, though there is some 
variation. Figure 14 shows the optimum buffer-velocity ratio as a function of the 
average-initial-velocity ratio Ug q. For higher initial-velocity ratios, the optimum 
buffer-velocity ratio is lower. In general, the optimum buffer velocities range from 
0. 2 to 0. 4 of the outer velocity for turbulent flow and from 0. 1 to 0. 3 for laminar flow. 
Thicker buffer regions require lower buffer velocities at moderate to high initial- 
velocity ratios. 

Figure 15 shows the effect of the thickness of the buffer region on containment 
factors that have been maximized with respect to the buffer-velocity ratio. As with 
buffer-velocity ratio, an optimum exists. Higher initial-average-velocity ratios re- 
quire somewhat thicker buffer regions, but the effect is relatively weak. For the range 
of conditions studied, the optimum momentum buffer thickness is approximately 1 jet 
radius. It is interesting to note that the optimum buffer thickness is virtually independent 
of whether the flow is laminar or turbulent. 
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CONCLUSIONS 


The results of this study indicate some characteristics of two- component, three- 
region coaxial flow. All the calculations were carried out for a channel- to jet-radius 
ratio of 4 and for the following physical properties: a ratio of inner- to outer-fluid 
molecular weight of 100, an inner-fluid Schmidt number of 1, a ratio of outer- to inner- 
fluid viscosity of 0. 2, and a ratio of outer- to inner-fluid molecular volume of 0. 2. For 
laminar-flow calculations, the jet Reynolds number was taken as 200; for turbulent flow, 
the jet Reynolds number was taken to be 20 000. Channel lengths of 4 and 8 jet radii were 
considered, and buffer thicknesses of 0, 0. 5, and 1 jet radius were investigated. Average 
ratios of outer- to inner-fluid initial velocity of 10, 50, and 100 were studied; for each 
of these cases, buffer- to outer-initial-velocity ratios of 0. 1, 0. 2, 0. 3, and 0. 4 were 
considered. For these ranges of conditions, the following results were obtained: 

1. An optimum buffer-region velocity exists that minimizes the acceleration of the 
inner stream by the faster moving outer stream; this maximizes the amount of inner 
stream gas contained within a given length of the outer duct. 

2. Similarly, there is an optimum buffer-region thickness. 

3. The optimum buffer velocities range from 0. 2 to 0. 4 of the outer velocity for 
turbulent flow and from 0. 1 to 0. 3 for laminar flow. 

4. The effect of a momentum buffer region is greater for higher ratios of outer- to 
inner-stream initial fluid velocity. 

5. For a duct radius that is four times that of the jet, the optimum buffer thickness 
is 1 jet radius for both laminar and turbulent flow. 

6. The presence of an intermediate-velocity momentum-buffer region of outer fluid 
significantly affects mass and momentum transfer and can increase the amount of inner 
fluid contained in a duct 4 jet radii in length by more than a factor of 2. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 17, 1965. 
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